Cosmic string parameter constraints and model analysis using small scale Cosmic 

Microwave Background data 
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We present a significant update of the constraints on the Abelian Higgs cosmic string tension 
by cosmic microwave background (CMB) data, enabled both by the use of new high-resolution 
CMB data from suborbital experiments as well as the latest results of the WMAP satellite, and by 
improved predictions for the impact of Abelian Higgs cosmic strings on the CMB power spectra. 
The new cosmic string spectra [1] were improved especially for small angular scales, through the use 
of larger Abelian Higgs string simulations and careful extrapolation. If Abelian Higgs strings are 
present then we find improved bounds on their contribution to the CMB anisotropics, /iq H < 0.095, 
and on their tension, G/xah < 0.57 x 10 -6 , both at 95% confidence level using WMAP7 data; and 
< 0.048 and G^ah < 0.42 x 10~ 6 using all the CMB data. We also find that using all the CMB 
data, a scale invariant initial perturbation spectrum, n 3 = 1, is now disfavoured at 2.4<r even if 
strings are present. A Bayesian model selection analysis no longer indicates a preference for strings. 



I. INTRODUCTION 

Observations of the cosmic microwave background 
(CMB) radiation [2-5] and the large-scale distribution of 
galaxies [6] are consistent with the inflationary paradigm 
(see e.g. [7]) in which cosmic structure was seeded in 
the very early stages of the universe. Attempts at rec- 
onciling observations with a late-time generation of the 
perturbations have not proven very successful and ap- 
pear to require a primordial mechanism like inflation [8- 
10]. Attention is now focused on the search for possi- 
ble clues on the detailed physical processes that created 
the initial perturbations, through tests for departures 
from the standard cosmological model, which is charac- 
terised by Gaussian adiabatic density perturbations with 
a power-law power spectrum. An important clue would 
be the detection of phase transition remnants like topo- 
logical defects, especially strings [11-16]. A large class 
of physically-motivated inflation models predict the ex- 
istence of cosmic strings, including hybrid inflation [17— 
21], brane inflation in the context of string theory [22, 23] 
and models rooted in supersymmetry and Grand Unified 
Theory (GUT) [24]. Previous analyses have shown that 
cosmological datasets leave room for significant effects 
from cosmic strings [25-28]. 

Data from the WMAP satellite, when interpreted via 
the standard ACDM model, suggests that the power 
spectrum of the primordial density perturbations obeys 
a power-law with index n s — 0.963 ± 0.015 [29]. This 
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appears to put some pressure on hybrid inflation models 
for which n s ~ 0.98, and is significantly different from 
the exactly scale invariant Harrison- Zel'dovich spectrum 
(n s = 1). However, previous work [25, 26] demonstrated 
that when cosmic strings are included the required value 
of n s can rise significantly, with the data then available 
implying n s — 1.00 ±0.03 with Abelian Higgs strings [25]. 
So long as an inflation model can provide the appropriate 
relationship between n s and the string tension, the data 
can instead be seen to favour the model, with examples 
being the D3/D7 brane inflation model [30] and semi- 
shifted hybrid inflation in the extended supersymmet- 
ric Pati-Salam model [31]. Additionally, if other types 
of cosmic defect such as semilocal strings [32] or cosmic 
textures [32, 33] are considered, then even larger n s val- 
ues are possible, with those models additionally offering 
lower values of the string tension or symmetry-breaking 
scale [32]. 

Here we revisit the constraints upon Abelian Higgs 
strings, incorporating updates both on the theory and 
the data side: We use more recent data, especially new 
observations on small angular scales (£ > 2000) where the 
string contribution has a greater relative importance [34] . 
The use of the small-scale data is enabled by a significant 
improvement in the CMB power spectrum predictions 
from Abelian Higgs string simulations [1], which bene- 
fits from improvements in computing power to achieve 
greater accuracy, and a better understanding of the gen- 
eration of perturbations to extend the spectra to a mul- 
tipole of I = 4000. 

There has been other recent work on cosmic string 
constraints [5, 35-37], with different string models and 
datasets. In [35] CMB data was taken from WMAP5 [29], 
ACBAR [3], BOOMERANG [38], CBI [39] and QUAD 
[40] experiments, and combined with galaxy power spec- 
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trum data from the SDSS Luminous Red Galaxy sample 
[41]. The data was fitted with the Unconnected Segment 
Model (USM) of cosmic string power spectra [42 -44], 
with parameters chosen to resemble Nambu-Goto strings 
(USM-NG) or an earlier Abelian Higgs power spectrum 
(USM- AH) [45]. The baryon density was constrained 
with a prior derived from estimates of the primordial 
Deuterium abundance in damped Lyman-a systems [46] . 
In [5] CMB data was taken from WMAP7 [2] and ACT 
[47], and fitted with the USM-NG cosmic string model, 
with priors on the Hubble parameter and BAO distance 
scales. The other parameters of the USM model beside 
the string tension can also be fitted to the data, rather 
than chosen to resemble numerical simulations [37] . The 
result is to weaken the constraints on the string tension. 
The apparent discrepancies between these works (includ- 
ing the present one) come both from the various datasets 
used and from the different input templates used for the 
string spectrum. The AH and the USM-NG cosmic string 
power spectra differ in shape, and more importantly in 
normalisation, leading to constraints on the cosmic string 
tension which are a factor of about 2 stronger for the the 
USM-NG model when fitted to the same datasets [35]. 

In this paper, we compare our latest Abelian Higgs 
string CMB anisotropy power spectra, exhibited in 
Fig. 16 of [1], to CMB temperature and polarization mea- 
surements, using 7-year WMAP data [2] augmented by 
small-scale part-sky data from ACBAR [3] and QUAD 
[4], and ACT [5]. Due to the possibility that the pres- 
ence of cosmic strings could introduce a scale dependent 
bias and induce early non-linear evolution in the mat- 
ter power spectrum, we prefer to be conservative and do 
not use the SDSS galaxy power spectrum data (see also 
section III). 

We present new constraints for the Abelian Higgs 
model on the parameter G/z, which is the product of New- 
ton's constant G and the string energy per unit length 
/j,, and the related parameter /io, which is the fractional 
power in the temperature power spectrum due to cos- 
mic strings at I = 10. We also discuss the implication 
of string presence upon the index n s , re-assessing the 
preference that we found for a 6-parameter cosmolog- 
ical model with Harrison-Zel'dovich perturbations and 
strings in [25]. We further consider the effect of adding 
data from measurements of the Hubble parameter (Hq) 
[48], baryon acoustic oscillations (BAO) [6] and the pri- 
mordial deuterium abundance (BBN) [46] as priors. We 
show that using CMB data alone, quite strong constraints 
can be derived without the uncertainty stemming from 
the use of other datasets, where the connection to the 
primordial perturbations is less easy to calculate. 



II. DATA AND METHODOLOGY 

Following the methodology used in [33] we add the 
defect CMB angular power spectra and the primordial 
CMB anisotropy angular power spectra (both tempera- 



ture and polarisation spectra) , and then explore the like- 
lihood of these combined Ci with the help of the Markov 
chain Monte-Carlo (MCMC) method, using a modified 
version of CosmoMC [49]. While we re-calculate the in- 
flationary contribution G™ f (h, n s , ...) at each point in pa- 
rameter space, we calculate the string contribution G| tr 
for a single set of cosmological parameters only. The 
normalization of string contribution to power spectrum 
is proportional to (G/i) 2 , where /i is the string mass per 
unit length and G is Newton's constant, so that 

C e (h,n s ,Gv,...) (1) 

= Cf{h,n s ,...)+ Cf(fto,G W ,..), 

where o denotes the values of the cosmological param- 
eters at which the string contribution was calculated: 
h = 0.72, n h h 2 = 0.214, and tt A = 0.75, and an opti- 
cal depth to last-scattering of r = 0.1. This approximate 
treatment of the string component significantly reduces 
the computational cost of the MCMC procedure and is 
justified as the allowed changes in cosmology affect the 
string component at a level below the observational un- 
certainties. 

We consider models both with and without strings and 
compare them to a number of datasets. To ease dis- 
cussing them, we shall denote the standard 6-parameter 
cosmological model, with the power spectrum of primor- 
dial curvature perturbations obeying a power-law in k, as 
model PL. The parameters are: the index of this power- 
law n s , its amplitude at k = 0.05 Mpc -1 A s , the Hubble 
parameter 100ft- kms~ 1 Mpc~ 1 , the physical baryonic and 
cold dark matter density parameters fl^h and fl c h 2 , the 
optical depth to last scattering r. As standard in Cos- 
moMC, we replace h with 9 which is approximately 100 
times the ratio of the sound horizon at last scattering to 
the angular diameter distance. As observations extend 
to smaller angular scales it is now common to include a 
seventh parameter Asz, the amplitude of the Komatsu 
and Seljak [50] template for the Sunyaev-Zel'dovich (SZ) 
effect, and we denote this model PLsz- Throughout our 
analysis we assume flat space, as would be yielded to a 
very good approximation by the relevant inflation mod- 
els. Special cases of PL and PLsz are obtained by set- 
ting n s = 1, the Harrison-Zel'dovich spectrum, giving 
"nested" models HZ and HZsz with 5 and 6 parameters 
respectively. 

To models PLgz and HZsz we add the contributions 
to the CMB power temperature and polarization spectra 
from Abelian Higgs cosmic strings (AH) with free nor- 
malization (G/iah) 2 , which we label as PLsz+AH and 
HZsz + AH respectively. We then compare the above 
mentioned models to seven-year WMAP data [2] (£ < 
1000), and study the effect of adding ACBAR [3] (300 < 
I < 3000) and QUAD data [4] (2000 < t < 3000), and 
finally ACT data [5] {£ < 10000). 

In order to fit ACT data, we also need to account 
for foreground point sources, which we model as Poisson 
noise with a normalisation A p [5]. As the point sources 
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FIG. 1. CMB data and predictions for the temperature angu- 
lar power spectrum from a model with strings. Parameters are 
those from the final column in Table I. The thick black line 
corresponds to the total power spectrum, the upper dashed 
black line corresponds to the inflationary component, the 
lower dashed black line corresponds to the string compo- 
nent, the humped grey line corresponds to the SZ contribution 
(from the ACT 148 GHz template) and the power-law grey 
line corresponds to the Poisson point source model. ACT 148 
GHz and ACBAR data are represented by red and blue bars 
respectively. The small glitch at t = 4000 is caused by the 
zeroing of the inflation and string components beyond this 
value. The result will be a slight overestimation of A p and 
therefore potentially a very slight underestimation of strings. 



are the dominant contribution for £ > 4000 we cut off 
the inflation and string power spectra for I > 4000. We 
show in Fig. 1 the data sets and the C/s used. 

We then investigate the effect of combining the CMB 
data with constraints from other measurements, which 
we impose as priors on the MCMC. We choose three of 
the most commonly used measurements: the Hubble pa- 
rameter (Hq) [48], the distance scales set by measure- 
ments of the Baryon Acoustic Oscillation (BAO) [6] and 
the baryon density inferred via standard Big Bang Nu- 
cleosynthesis from measurements of the primordial deu- 
terium abundance (BBN) [46]. 



III. PARAMETER CONSTRAINTS 

A. CMB data only 

We initially consider the comparison of the models 
PLsz and PLgz+AH against (only) CMB data, with re- 
sults shown in Table F Previously, when we performed a 
similar analysis to the three- year WMAP data we found a 
preference for a non-zero contribution from cosmic strings 
and an almost exactly scale-invariant power spectrum: 
n s = 1.00 ± 0.03 [25]. We found that with WMAP7 the 



tendency remains, and both n s and Sl^h 2 tend to higher 
values when a string contribution is added. 

This increase in n s upon the inclusion of cosmic de- 
fects is due to a parameter degeneracy noted in Ref. [33] . 
When a contribution from cosmic defects is added to the 
adiabatic scalar component seeded by inflation, it is clear 
that A s must be reduced if the fit to the WMAP data is 
to be maintained. However, the inflationary contribution 
is most sensitive to A s at the first peak and therefore this 
peak is excessively lowered. This can be countered by an 
increase in flhh 2 , which lifts the first peak while lower- 
ing the second. Additionally n s can be raised to lift the 
second peak relative to the first. Also, an increase in h 
helps reverse the shift of the first peak to smaller scales 
which is caused by the latter, and also the shape of the 
string component. 

However, the increase in n s is not as high as with 
WMAP3, and the exactly scale-invariant power spectrum 
is not as favoured with the newer WMAP7. The mean 
n s value is reduced to n s = 0.985 ± 0.020, and therefore 
there is now a slight preference for n s < 1, even when 
local cosmic strings are included. Hence inflation models 
that yield n s values greater than 1 cannot be as easily 
aided by the inclusion of cosmic strings. 

If we now add ACBAR and QUAD at high £, then 
we find additional constraints upon the string compo- 
nent, reducing its mean amplitude, as well as a lower 
mean value for n s . Finally, further adding ACT data, 
the constraints on the string component get even more 
stringent, and the value of n s is now hardly shifted up- 
wards when we allow for a string contribution. Use of the 
ACT data however necessitates the addition of SZ and 
point-source high-f templates. Here we only quote re- 
sults when allowing for the presence of unclustered point 
sources, which give rise to an extra contribution with con- 
stant Cg . This population is detected at high significance 
(A p ss 17 ± 2, see tables I and II). We also experimented 
with adding a clustered point-source population [5] , find- 
ing that the amplitude of the template power spectrum 
was A c « 3 ± 3. As the result is consistent with zero, 
we did not include clustered point sources in our main 
results. 

Table III shows a comparison of 95% upper bounds on 
the string tension obtained by different data sets and dif- 
ferent ways of modelling the string spectra. To date, the 
power spectra have been calculated using classical field 
theory simulations of the Abelian Higgs model [1, 45] or 
by the Unconnected Segment Model (USM), which repre- 
sents the evolving string network with moving sticks with 
mass and tension. The sticks have random positions and 
velocities, and disappear so that their overall energy den- 
sity scales as i~ 2 . There is enough freedom in the USM 
to model the important scalar and vector components 
of the Abelian Higgs power spectra (USM- AH), and to 
make a prediction for the power spectra Nambu-Goto 
simulations would produce based on their mean square 
velocity and density (USM-NG). Note that there is cur- 
rently no calculation of the CMB power spectra using NG 
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Data 


WMAP7 


+ACBAR+QUAD 


+ACT 


Model 


PL S z 


PL SZ + AH 


PLsz 


PL SZ + AH 


PLsz 


PL SZ + AH 


Q h h z 


0.0224(6) 


0.0235(10) 


0.0225(5) 


0.0231(6) 


0.0223(5) 


0.0226(5) 


n c h 2 


0.111(5) 


0.101(6) 


0.112(5) 


0.109(5) 


0.113(5) 


0.111(5) 


9 


1.039(3) 


1.042(3) 


1.041(2) 


1.042(2) 


1.041(2) 


1.042(2) 


T 


0.085(15) 


0.094(17) 


0.088(15) 


0.086(15) 


0.087(14) 


0.085(13) 


n B 


0.967(13) 


0.985(20) 


0.969(13) 


0.975(13) 


0.965(12) 


0.969(13) 


ln(10 10 A s ) 


3.07(4) 


3.07(5) 


3.08(3) 


3.05(3) 


3.09(3) 


3.07(3) 


^sz 


1.1(6) 


0.9(6) 


1.1(5) 


0.9(6) 


0.5(3) 


0.4(3) 












16(2) 


17(2) 


10 12 (G/x) 2 




0.13(10) 




0.11(7) 




0.07(5) 


h 


0.708(23) 


0.735(38) 


0.712(23) 


0.731(25) 


0.703(23) 


0.717(25) 


fw 




0.035(30) 




0.031(20) 




0.020(14) 


Likelihood 


3734.044 


3733.997 


3807.694 


3807.321 


3823.176 


3823.018 



TABLE I. Marginalized likelihood constraints on model parameters. The basis parameters, upon which we assume flat priors, 
are listed above the horizontal line, with other derived parameters listed below this line. The notation is such that 0.0224(6) = 
0.0224 ± 0.0006, i.e., it indicates an expectation value of 0.0224 and a standard deviation uncertainty of 6 in the final digit. 
In the first two columns, the data used are those from WMAP7. The next two columns use ACBAR and QUAD, as well as 
WMAP7. The last two columns use all the data from WMAP7, ACBAR, QUAD and ACT. Since these last two columns use 
ACT, we also include a A p to take into account Poisson point sources. 



Data 


H 


BAO 


BAO+flo 


BBN 


BBN+if<, 


BAO+BBN 


BAO+BBN+ifo 




0.0228(5) 


0.0224(4) 


0.0225(4) 


0.0226(5) 


0.0227(5) 


0.0224(5) 


0.0225(5) 


n c h 2 


0.110(3) 


0.114(3) 


0.113(3) 


0.111(5) 


0.110(4) 


0.114(3) 


0.113(3) 


e 


1.042(2) 


1.041(2) 


1.041(2) 


1.041(2) 


1.042(2) 


1.041(2) 


1.041(2) 


T 


0.087(14) 


0.083(13) 


0.085(14) 


0.088(15) 


0.090(15) 


0.085(14) 


0.086(14) 


n s 


0.972(11) 


0.964(11) 


0.966(10) 


0.968(12) 


0.972(12) 


0.963(11) 


0.966(11) 


ln(10 10 A s ) 


3.06(3) 


3.07(3) 


3.07(3) 


3.07(3) 


3.07(3) 


3.08(3) 


3.08(3) 


A sz 


0.4(3) 


0.4(3) 


0.90(6) 


0.4(3) 


0.4(3) 


0.4(3) 


0.4(3) 


A 


17(2) 


17(2) 


17(2) 


17(2) 


17(2) 


17(2) 


17(2) 


10 12 (G M ) 2 


0.08(5) 


0.06(4) 


0.07(5) 


0.07(5) 


0.08(5) 


0.06(4) 


0.06(5) 


h 


0.725(21) 


0.701(15) 


0.707(14) 


0.714(24) 


0.723(20) 


0.700(15) 


0.706(14) 


ho 


0.022(15) 


0.016(12) 


0.017(12) 


0.019(13) 


0.020(14) 


0.016(12) 


0.017(12) 



TABLE II. Marginalized likelihood constraints on model parameters, as in Table I. The notation is such that 0.0228(5) = 
0.0228 ± 0.0005, i.e., it indicates an expectation value of 0.0228 and a standard deviation uncertainty of 5 in the final digit. 
In this case we fit PL S z+AH with all the CMB data (WMAP7, ACBAR, QUAD, ACT), with different priors. The different 
priors include (combinations of) priors coming from determinations of the Hubble parameter (-ffo), Baryon Acoustic Oscillations 
(BAO) and Big Bang Nucleosynthesis (BBN), as explained in Section IIIB. 



simulations and the Unequal Time Correlator method of 
[!]■ 

In general, we see that the WMAP7 data puts tighter 
constraints on the string contribution, mostly due to its 
smaller errors on the data points. The extra high-£ CMB 
experiments further push down the string contributions, 
where strings are expected to be more important. How- 
ever, even though the constraints are tighter, B-mode ex- 
periments such as CMBpol will be able to detect string 
components that are much lower than the bounds ob- 
tained here [51]. 

Note that the USM gives tighter constraints when 
modelling Nambu-Goto strings. This is mainly due to 
two features of NG simulations: a higher string density 
because of the absence of the decay channel into massive 
radiation [52] , and the persistence of the higher radiation 
era string density well into the matter era [36] . 

Figure 2 shows marginalised 2-D likelihoods of some 
of the parameters in the models, and depicts how the 



constraints have shrunk by taking into account the new 
data. The contours obtained with WMAP3 data were 
much wider, and it can readily be seen that n s = 1 was 
close to the mean value. The degeneracies among the 
parameters are also apparent. The inclusion of WMAP7 
data made the contours smaller by more than a factor of 
two, leaving less space for strings, and pushing n s away 
from scale invariance. By including the rest of the CMB 
data mentioned above, the contours shrink an extra fac- 
tor of two. 



B. Adding non-CMB data 

In this section we investigate the effect of combin- 
ing the CMB data with constraints from other mea- 
surements, the Hubble parameter (H ) [48], the distance 
scales set by measurements of the Baryon Acoustic Os- 
cillation (BAO) [6] and the baryon density inferred via 
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Model 


Data set 


10 b Gfi (95%) 


fio (95%) 


AH [25] 


WMAP3+BOOMERANG+CBI+ACBAR+VSA 


0.7 


0.11 


AH (this work) 


WMAP7 


0.57 


0.095 


AH (this work) 


WMAP7 + ACBAR + QUAD + ACT 


0.42 


0.048 


USM-AH [35] 


WMAP5 


0.68 


0.11 


USM-NG [35] 


WMAP5 


0.28 


0.054 


USM-NG [5] 


WMAP7+ACT 


0.16 





TABLE III. CMB constraints on string tension obtained by different string modeling, and different data sets. The table shows 
95% upper bounds on G/i. The different models shown are the Abelian Higgs model (AH) and the Unconnected String Segment 
model (USM), with parameters chosen to resemble the Nambu-Goto (NG) or Abelian Higgs (AH) spectra. The difference in 
the values of G\x and /io are due to both the different string templates used and the different data sets used. 




FIG. 2. Marginalised 2-D likelihood contours for different 
data sets. The unshaded regions in both the top and bot- 
tom panel correspond to the 1-a and 2-a contours obtained 
in [25] when contrasted with WMAP3 data [53]. The shaded 
regions correspond to the 1-a and 2-a contours obtained in 
the present work. The top panel shows the contours for 
WMAP7 data, and the bottom panel shows the contours for 
WMPA7+AQBAR+QUAD+ACT. 



standard Big Bang Nucleosynthesis from measurements 
of the primordial deuterium abundance (BBN) [46] . We 
include the data as priors on the MCMC. 

We see that the most recent determination of Ho of 
74.2 ±3.6 kms -1 Mpc -1 from a combination of Cepheids 
and SN-Ia data [48] is slightly higher than the CMB mean 
value of 70.3 ± 2.3 kms -1 Mpc -1 under the PLsz model. 
Imposing this prior on H for the PLsz+AH model al- 
lows for a slightly higher value of /io, but the effect is 
negligible. 

Distance measurements with the help of the baryon 
acoustic oscillations (BAO) are supposed to be very re- 
liable as they are based on a physical mechanism closely 
related to the one that created the acoustic oscillations 
in the CMB. However, the resulting oscillatory feature 
in the galaxy power spectrum is on intermediate scales 
where the matter power spectrum P(k) is on the verge of 
deviating from the linear prediction already in the stan- 
dard cosmological model without strings. In addition, 



measurements of the galaxy power spectrum or corre- 
lation function are always subject to galaxy bias. The 
reliability of the BAO measurement has been tested with 
large N-body simulations in the absence of topological 
defects [54] but if cosmic strings are present then one 
would need to check whether they are able to affect ei- 
ther the galaxy bias or the scale on which perturbations 
turn non-linear. This requires the inclusion of string per- 
turbations in N-body simulations [55]. In any case the 
effect of adding the BAO data is to shift mean values by 
less than half a standard deviation. 

The noted increase in tl^h 2 for PL+AH is of inter- 
est since model PL already yields values of this parame- 
ter which are larger than those measured by alternative 
means. Big bang nucleosynthesis (BBN) calculations link 
the primordial abundance of the light isotopes to flbh 2 
but it is only possible to measure such abundajnees in 
astronomical objects at more recent times. While ilb/i 2 
measures using different isotopes do not wholly agree, 
deuterium abundances in high redshift gas clouds (z « 3) 
can be measured via absorption lines in quasar spectra, 
yielding a result of Q h h 2 = 0.0213 ± 0.0010 [46]. Assum- 
ing a Gaussian uncertainty, this can be included into our 
analysis, yielding the results in Table II. However, given 
that there are only 7 measurements, with abundances 
corresponding to Vl^ti 2 in the range 0.017 to 0.030, the 
BBN constraint should be treated with caution. The ef- 
fect on model PLsz is a small reduction (less than half a 
standard deviation) in Q,\,h 2 via a degeneracy similar to 
that noted for the defect case above. The effect on the 
parameters n s and /io of the model PLsz+AH is negli- 
gible, showing that CMB data are now good enough to 
render the BBN constraint unnecessary. 

For these reasons we still consider the CMB-only con- 
straints as the most reliable results, which we therefore 
use as the main figures quoted in this paper. In any case, 
the addition of the small-scale CMB data is as powerful 
as the addition of the non-CMB data considered here. 



IV. MODEL SELECTION 

We now compare models via maximum likelihood val- 
ues £ ma x and also via Bayesian evidence ratios, which 
yields a more complete analysis of the freedom within a 
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Data 


CMB only 


Measure 


Axcff 


InB 


PL SZ 








PLsz + AH 


-0.3 


-3.3 + 0.1 


HZsz 


+7.8 


-1.0 + 0.3 


HZsz+AH 


+4.7 


-3.2 + 0.2 



TABLE IV. A\ 2 and the natural logarithm of the Bayes fac- 
tors (In B) for four different models, taking as the base model 
PLsz- The data used in these comparisons are the full CMB 
data used in this work, which means that ACT is also taken 
into account, and the models have foregrounds modelled with 
Poisson point sources. Therefore, PLsz is an 8 parameter 
model, with the usual 6 parameters plus Asz and A p . 

model (see, for example, [56]). We determine maximum 
likelihood values via low temperature MCMC chains and 
then express results via Ax 2 ff = — 2Aln(£ max ) between 
two models. We determine evidence ratios using the 
Savage-Dickey density ratio [57, 58], which can be used 
with nested models. This yields the Bayes factor between 
an n+1 parameter model M n +i (such as PLsz+AH) and 
an n parameter subset of that model M. n (such as PLsz), 
via: 

_ P(d\M n+1 ) _ ir(x = x*\d,M n+1 ) 

*n+v»- P{d \ Mn) P (x = x*\d,M n+i y ( ) 

where x is the (n + l)th parameter and x* its fixed 
value in the n parameter model. That is, the complex 
calculation of evidences is reduced to determining the 
value of the (independent) prior upon x in model M. n +\, 
it(x = \d, M n +i), and the marginalized posterior prob- 
ability density for x at x* (given the data d and model 
M n+ i) P(x = x*|d, M. n+ \). The former is trivial for the 
flat priors that we assume here but highlights the depen- 
dence of the results upon prior choice. The latter, on the 
other hand, can be found via binning the MCMC results 
according to their x values. 

Our results for both statistics are shown in Table IV for 
the models employed above, based around adding an AH 
component to model PLsz- We now also include results 
for model HZsz (with and without an AH component). 
These models are nested into each other at points n s = 1 
and (Gfj,) 2 — 0. There are hence two flat priors whose 
range is relevant to our results and for these we take: 
n s = 0.75 -> 1.25 and (Gfi) 2 = -> (2 x 1CT 6 ) 2 . The 
latter is inspired by the approximate limit that we would 
obtain if we compared our string spectra to COBE data 
(£< 10) l . 

While in [32] we found that the addition of an AH 
string component to the PL model lowered the x 2 by 
3.5 and lead to a comparable model probability for 



This is roughly equivalent near the nesting point /io = to the 
prior of /io = — > 1 that we employed in Ref . [25] , although then 
the limits were simply mathematical constraints on a fractional 
quantity and near /io = l the two priors would noticeably differ. 



0.07 
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/ 

FIG. 3. Cumulative model-averaged ID probability distri- 
bution for /io using the Bayes factors in table IV. For the 
models and priors chosen, 94.6% of the weight is in the "no 
strings" models PLsz and HZsz. From the choice of priors 
and models, this figure gives an idea of the constraints on 
generic scenarios where Abelian-Higgs strings are formed at 
an energy scale associated with inflation (see text for more 
details on priors and model choices). 



PL+AH (and even a somewhat higher model probabil- 
ity for HZ+AH) we now find that all alternative models 
are weakly to moderately disfavoured when compared to 
PLsz- There are two reasons for the change: Firstly 
the improvement of fit when adding AH to PLsz has de- 
creased, with HZsz+AH now disfavoured, and secondly 
the posterior parameter volume (i.e. the error bar on 
(Gfj,) 2 ) has also decreased, which increases the "Occams 
razor" penalty for a fixed prior size. The change in good- 
ness of fit (due to n s = 1 being disfavoured now even 
when strings are present) is especially important for the 
HZsz+AH model. 

Given the relative model probabilities, we are able to 
take into account the model uncertainty in the determi- 
nation of (Gfj,) 2 or /io by marginalising over the models 
(see e.g. [59]): 

P(fio\d) = J2 p (ho,M\d) = J2P(fw\d,M)P(M\d). 

M M 

The sum over models here runs over the set {PLsz, HZsz, 
PLsz+AH, HZsz+AH }. If we take the prior probability 
of all models to be equal then we can replace the model 
probability P(A4\d) with P(d\A4), and since we deter- 
mine the overall normalisation of the model-averaged 
posterior by requiring that the integral over /io be equal 
to 1 we can in addition replace P(d\M) with Bayes fac- 
tors relative to one of the models. 

Since for the two models without strings we have that 
-P(/io) = <5(/io), which is difficult to plot, in Fig. 3 we 
show the cumulative probability that /io is larger than a 
given value. Based on the Bayes factors in table IV, we 
find that 94.6% of the posterior volume is in the two "no- 
string" models PLsz and HZsz- The use of the model- 
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averaged posterior would therefore tighten the 95% up- 
per limit on strings significantly. However, this result is 
strongly dependent on the selection of models considered, 
and we prefer to quote upper limits on /i and Gfi under 
the assumption that Abclian-Higgs strings are present 
and that the initial power spectrum is of a power-law 
form, i.e., the model PLsz+AH. 

Note that the Bayes factors in Table IV also de- 
pend strongly on the priors used, especially on the fact 
that we assume G[i to be uniformly distributed in the 
range [0,2 x 1CP 6 ]. This prior appears appropriate for 
"GUT strings" associated with the energy scale of infla- 
tion. If we are interested in constraints on strings that 
were formed in a phase transition at an arbitrary energy 
scale, then a flat prior in lnGjti would be better moti- 
vated. Such an alternative prior would change the result 
strongly, as the prior weight is then equally distributed 
over the range of possible scales for GpL (a Jeffreys prior 
[60] ) 2 . Such a prior, with a sensible choice of the lnG/z 
lower bound, would place a large prior weight on very low 
G[i values where the posterior is flat. As the posterior 
is close to maximal for f w — > with current data it is 
clear that a strong preference for or against strings is im- 
possible to obtain in both model selection and the above 
parameter results, resulting in roughly equal model prob- 
abilities for PLsz and PLsz+AH. We refer the reader to 
the detailed discussion of the computational methods and 
the expected effects of a logarithmic prior on G\i in [51]. 

V. CONCLUSIONS 

In this paper we use an improved determination of the 
Abelian Higgs cosmic string CMB power spectrum Ci up 
to t = 4000, combined with new CMB data that extends 
out to I w 10000 to update our constraints on the maxi- 
mal contribution of strings to the CMB anisotropies. We 
find that /A/ 1 < 0.048 at 95% confidence level, which 
translates into a bound on the Abelian Higgs string ten- 
sion of G/iAH < 0.42. Using only WMAP7 data, the 
bounds we obtain arc /j^ H < 0.095 and G^ah < 0.57 
at 95% confidence level. These bounds are somewhat 
weaker than the ones derived by other groups using the 
unconnected segment model for Nambu-Goto strings, 
which is mostly due to a lower string density in the 
field theory simulations. Note that even though the con- 
straints are tighter than previous ones, the Planck satel- 
lite is expected to improve them (to about /io < 0.01) 
[61], and a future dedicated B-mode satellite will be able 
to detect string components that are much lower than 
the bounds obtained here [37, 51]. 



While the WMAP 7-year data set already improves 
the constraints significantly compared to the 3-year data 
set, the additional constraints from high-^? data (ACBAR, 
QUAD and especially ACT) help to break the degeneracy 
between the amplitude, n s and Slbh 2 so that a Harrison- 
Zel'dovic spectrum with n s = 1 is now also disfavoured 
when allowing for the presence of strings. Because of this, 
we now obtain strong constraints with CMB data alone, 
and the addition of BAO, BBN and Hq constraints is less 
important than in the past. 

Previously Bayesian model selection slightly favoured 
HZ+AH over PL when using only CMB data [25, 32]. 
This has changed with the improved CMB data that no 
longer prefers a non-zero /j^, H and also puts pressure on a 
Harrison-Zel'dovich initial power spectrum even with AH 
strings. For the choice of priors and models investigated 
here we find that models with AH strings are disfavoured. 

We expect that the strong constraints obtained in this 
work put pressure on generic models that form strings 
of the Abclian-Higgs type at energy scales associated 
with inflation, e.g., these constraints push hybrid infla- 
tion models to small couplings. We postpone the study 
of the constraints on specific implementations of high en- 
ergy physics based models of inflation that form defects 
to later work. Generically, in those models, the string 
tension is not a free parameter, since it is related to other 
measurable predictions. Having to allow for a low string 
tension will put pressure on other parameters to meet the 
measurements, and it will be interesting to investigate 
whether the interplay between the parameters is able to 
render the model viable. 
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